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ABSTRACT 

We discuss the cosmological implications of the new constraints on the power spectrum of the Cosmic Mi- 
crowave Background Anisotropy derived from a new high resolution analysis of the MAXIMA-I measurement. 
The power spectrum indicates excess power at I ~ 860 over the average level of power at 41 1 < I < 785. This ex- 
cess is statistically significant on the ^95% confidence level. Its position coincides with that of the third acoustic 
peak as predicted by generic inflationary models, selected to fit the first acoustic peak as observed in the data. The 
height of the excess power agrees with the predictions of a family of inflationary models with cosmological param- 
eters that are fixed to fit the CMB data previously provided by BOOMERANG-LDB and MAXIMA-I experiments. 
Our results, therefore, lend support for inflationary models and more generally for the dominance of adiabatic 
coherent perturbations in the structure formation of the Universe. At the same time, they seem to disfavor a large 
variety of the non-standard (but inflation-based) models that have been proposed to improve the quality of fits to 
the CMB data and consistency with other cosmological observables. 

Within standard inflationary models, our results combined with the COBE-DMR data give best fit values and 
95% confidence limits for the baryon density, Vt b h 2 ~ 0.033±0.013, and the total density, n = 0.9 The 
primordial spectrum slope (n s ) and the optical depth to the last scattering surface (r c ) are found to be degenerate 
and to obey the relation n s ~ (0.99±0.14) + 0.46r c , for r c < 0.5 (all 95% c.l.). 

Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

Observations of Cosmic Microwave Background (CMB) 
temperature anisotropy are reaching maturity. The high signal- 
to-noise multi-frequency data gathered by the balloon-borne 
BOOMERANG-LDB (de Bernardis et al. 2000) and MAXIMA-I 
(Hanany et al. 2000) experiments set stringent constraints on 
the shape of the power spectrum in the broad range of angular 
scales ranging from ~ 5° down to ~ 10' scales (corresponding 
to a range in £-space from ~ 50 up to ~ 600 for BOOMERANG- 
LDB and from ~ 35 to ~ 800 for MAXIMA-I. The measure- 
ments firmly established the existence of a peak in the power 
spectrum at £ <~ 220, also suggested by the data from earlier 
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observations (e.g., Miller et al. 1999 and Mauskopf et al. 2000). 
Because no secondary peaks were indisputably seen, there was 
no unambiguous evidence for an inflation-like scenario where 
structure formation is driven by passive adiabatic coherent fluc- 
tuations (but see also results concurrent with ours - Netterfield 
et al. 2001, and Halverson et al. 2001). 

The level of power detected beyond the first peak, on subde- 
gree scales, was found to be somewhat lower than that generally 
expected for inflation-based "concordance" models of pre-2000 
(e.g., Ostriker & Steinhardt 1995, Krauss & Turner 1995). Al- 
though very good fits to the MAXIMA-I and BOOMERANG-LDB 
data could be found within the inflationary family of models 
(e.g., Jaffe et al. 2001), the most favored values of the baryon 
density were shown to be higher than that inferred from the 
standard arguments based on the cosmological (Big Bang) nu- 
cleosynthesis (BBN), though the latter was found to be within 
the 95% (97%) confidence limits derived from the MAXIMA-I 
(maxima-i+boomerang-ldb) measurement. Notwithstand- 
ing this fact a number of alternatives/extensions to the standard 
cosmologies were suggested (e.g., Bouchet et al. 2000, Enqvist, 
Kurki-Suonio, & Valiviita, 2000, Peebles, Seager, & Hu, 2001). 
Further high resolution data are required to test these models. 

A companion Letter, Lee et al. (2001), presents a new anal- 
ysis of the MAXIMA-I data that extends the measured power 
spectrum from the previously published range 35 < I < 785 up 
to £ < 1235. In this Letter we discuss the cosmological signifi- 
cance of the new result. 

This Letter is organized as follows. In section 2, we look for 
the signature of acoustic oscillations in our power spectrum, 
the existence of which is predicted by inflation-motivated sce- 
narios. In section 3, we discuss the cosmological implications 
of such a feature and derive constraints on some of the cos- 



1 



2 



Stompor et al. 



mological parameters within the family of standard inflationary 
cosmological models. 
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Fig. 1. — Angular power spectra of the CMB anisotropy recovered 
from the high resolution map of the MAXIMA-I (Lee et al. 2001). The 
shaded rectangles show the ranges of the bins used in our likelihood 
analysis in Section 2, and the 68% confidence limits on the power level 
in these bins. The dashed line shows the best-fit model to the previous 
low resolution results of MAXIMA-I (Hanany et al. 2000), the solid 
line shows the best-fit model for the new data (Lee et al. 2001), and 
the dotted line shows the best-fit model, fulfilling the BBN constraint, 
found within the subfamily of all models considered here. 

Finally, in section 4 we present our conclusions and comment 
briefly on the viability of non-standard cosmological models in 
light of the new data. 

2. THE EXCESS POWER AT £ ~ 860 ± 75 

Lee et al. (2001) set useful constraints on the CMB 
anisotropy power spectrum up to t ~ 1200, corresponding to 
angular scales down to 5'. In terms of the simple x 2 statis- 
tics 1 , this power, given as 1(1 + l)Q/27r in Figure 1, is consis- 
tent with a straight horizontal line for I £ 410 yielding a value 
X 2 — 5.2 for 7 degrees of freedom. However, not only does 
such a model have no physical justification, but also the data 
themselves qualitatively suggest the presence of high power at 
i ~ 860, in excess of the power level at I ~ 410-785, and 
£ ~ 1000- 1200. In fact a feature at I ~ 850, a third acoustic 
peak, is anticipated within a family of generic inflation-based 
models which are selected to reproduce a first peak at a position 
I ~ 220 as seen in the data. A statistically significant detection 
of such a feature would provide important additional support 
to this family of models. Statistical and systematic errors and 
correlations between the spectral bands obscure the statistical 
significance of the results. We therefore employ a likelihood 
approach to assess the statistical importance of the feature at 
I ~ 860 in the MAXIMA- 1 power spectrum. 

We work within a three-dimensional parameter space. The 
initial parameters are the bin powers Co, C\ and C2, for bins, Bo, 
411 < I < 785, B u 786 < I < 925, and, B 2 , 926<£< 1235, as 
shown by the shaded regions in Figure 1 . Those are chosen to 
facilitate the search for the feature at I ~ 850. Such a family of 
spectrum models includes a flat spectrum as a special case, and 
in the Bayesian spirit our analysis seeks to determine how well 
such a model represents our data. 

'Here and below, we define \ 2 assuming a Gaussian distribution rather than 




Fig. 2. — Left: Two-dimensional likelihood as a function of the bin 
power ratios. The contours correspond to 68, 95 & 99% confidence 
levels computed through integration over the regions enclosed within 
given contours. The shading shows the area, which correspond to the 
existence of the "excess power" in the bin centered at I ~ 860. Right: 
The one-dimensional likelihoods for bin power ratios TZm. Different 
lines are used for the two definitions of Bq. The solid lines are for 
41 1 < £ < 785 and the dashed lines for 486 < £ < 785. 

The likelihood for the three-dimensional parameter space 
is computed assuming the offset lognormal approximation of 
Bond, Jaffe & Knox (2000) to the probability distribution of 
the bin powers, denoted hereafter L (Co, C\ ,€2)- The likelihood 
also depends on "nuisance" parameters describing the contri- 
bution of systematic effects such as the total calibration uncer- 
tainty and the beam and pointing error. We model all of these 
as fully correlated between bins and Gaussian-distributed with 
a dispersion depending on I as given in Lee et al. (2001). We 
marginalize over these parameters by numerical integration. 

The question we first ask is whether the power in the two 
side bins (Bq and B2) is lower than the power in the middle bin 
centered at I <~ 860, and at what confidence level. For that rea- 
son, we introduce two new parameters, 72oi and 7v-2i, which are 
given by the ratio of power between bins IZij = Cij Cj . As a third 
parameter, describing the overall normalization, we choose the 
bin power amplitude defined above: C\ . We adopt flat priors in 
the bin powers and marginalize over the parameter, C\, leaving 
only a two dimensional problem, 

£(Ro U K 2i ) oc J dd C 1 2 £[C (7e i,Ci),Ci,C 2 (7e 2 i,Ci)], 

(1) 

where C\ is a Jacobian. The results are shown in Figure 2. The 
left panel shows the contours of the two dimensional likelihood 
£(7^.01,7^21) as a function of the ratios 7?-oi and 72-21- The con- 
tours show 68, 95, 99% confidence levels computed assuming 
flat priors for both parameters. The shaded area indicates the 
region of parameter space favoring excess power in the middle 
bin relative to its neighbors. Models without excess power can 
not be rejected at a confidence level higher than "la" (~ 68%). 

It is also apparent, however, that most of the models pre- 
ferred by the data, have power in the leftmost bin, Bo, that is 
lower than that in the central bin, B x . This point is addressed in 
the right panel of Figure 2, showing the one-dimensional like- 
lihoods computed through an explicit marginalization over the 
other parameter. We find that 1Zq\ = 0.49^'^ and hence it is 
lower than 1 at the confidence level ~ 95% (cf., the value of 
72-01 computed in the next section), yielding a "2<r" detection 
of the power rise at £ ~ 860 over the intermediate-^-range cov- 
ered by the Bo bin. This result remains unchanged if we allow 

approximately offset-lognormal form used in the explicit likelihood calculations. 



Cosmological implications of the MAXIMA-I CMB anisotropy measurement. 
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for the presence of the (uncorrelated) point-source-like compo- 
nent with the amplitude as expected for the 150 GHz band of 
MAXIMA-I even in the most pessimistic cases (Lee et al. 2001). 

It may appear that this high confidence level is largely due 
to the leftmost point of the spectrum (at I ~ 435) included in 
our (clearly arbitrary) definition of Bq. However, excluding that 
point lowers the confidence level only to ~ 93%. Because of 
the large uncertainties of the Lee et al. results beyond I ~ 900 
the statistical confidence that the power declines beyond the bin 
B\ is only ~ 80%. However, the recent CBI result (Padin et al. 
2001) constrains the power to 882^ [iK 2 at I ~ 1 190^ , pro- 
viding extra and independent support to the presence of such a 
decline. The high power level we see in the MAXIMA- 1 data at 
I ~ 860, therefore can not extend far beyond the right edge of 
the bin B\ (~ 925). The amplitude of this excess power is re- 
stricted by our data to C\ = 3273 l^tiK 2 at 95% c.l. including 
both statistical and systematic errors. 

In the usual family of inflation-based models, one also ex- 
pects to find a second peak in the region of 410 < t < 785. 
While a single flat band power provides an excellent fit to all 
the points within bin Bo with X 2 ~ 1-5 and 0.26 for 4 (wider 
bin) and 3 (narrower) degrees of freedom 2 , the presence of a 
second peak in that region is still comfortably admissible by 
the data (Figure 1). With our choice of the range of Bo, this bin 
is expected to be only a factor £ 2 wider than a typical feature 
of the power spectrum (A£ ~ 150) of the standard inflationary 
model. If such a peak structure is present in that range, /?oi < 1 
means that the bin power at I ~ 860 is higher than the level of 
power in some subsection of the 410 < I < 786 range. Then 
our analysis above demonstrates that this happens at a confi- 
dence level of at least 95% and that value is our confidence 
level of detecting a feature in the power spectrum at I ~ 860. 

3. COSMOLOGICAL PARAMETERS 
3.1. Parameter space 

Figure 1 demonstrates that an inflationary model with pa- 
rameters as determined by fitting to the MAXIMA-I power spec- 
trum published by Hanany et al. (2000) provides a very good 
fit over the entire current range of our new data together with 
the COBE-DMR result. The total \ 2 is ^ 32 for 41 data 
points (MAXIMA-I plus DMR) and \ 2 ~ 8 for the 13 points 
of MAXIMA-I only. Thus, our new data are consistent with 
the predictions made on the basis of the relatively narrow class 
of inflation-motivated models assumed in the previous papers. 
Such a class of models can be seen as a particular subset of the 
family of models considered in the previous Section. 

Here we consider a seven-dimensional space of parameters. 
The parameters include the amplitude of fluctuations at I = 10, 
Cio, the physical baryon density, fl\,h 2 , the physical density of 
cold dark matter, Vt c & m h 2 , and the cosmological constant, f^A, 
the total energy density of the universe, ft = f^b + ^cdm + ^A, 
the spectral index of primordial scalar fluctuations, n s , and the 
optical depth of reionization, r c . We use the following ranges 
and sampling: Cm is continuous; fl = 0.3, 0.5, 0.6, 0.7, 0.75, 
1.2, 1.3, 1.5; fl b h 2 = 0.00325, 0.00625, 0.01, 0.015, 0.02, 
0.0225, 0.04, 0.045, 0.05, 0.075, 0.1; n cdm h 2 = 0.03, 0.06, 
0.12, 0.17, 0.22, 0.27, 0.33, 0.40, 0.55, 0.8; Q A = 0.0, 0.1, 0.2, 
1.0; « s = 0.6, 0.7, 0.75, 0.8, 0.85, 0.875, 1.2, 1.25, 1.5; 
r c = 0, 0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.3, 0.5. The justifica- 
tion of the presented choice of the parameter space can be found 

2 Due to the strongly asymmetric shape of the \ 2 distribution with so few 
distributions. 



in Balbi et al. 2000. Models were computed using a version of 
CMBFAST by Tegmark, M., Zaldarriaga, M., & Hamilton, A., 
(2001), originally by Seljak & Zaldarriaga (1996). 

3.2. Results 

We compute the likelihood on the grid for models using an 
offset lognormal approximation (Bond, Jaffe, & Knox 2000), 
including both statistical and systematic errors in a manner 
analogous to Sect. 2. We neglect a subdominant pointing un- 
certainty (Lee et al. 2001). The likelihood for a subset of pa- 
rameters is evaluated by an explicit marginalization over all the 
remaining parameters, using a top-hat prior in parameters used 
to define our database. In addition, we impose top-hat priors 
for the value of the Hubble constant 0.4 < h < 0.9, the mat- 
ter density, fl m = fib + ^cdm > 0.1 and the age of the Universe 
> 10 Gyr. We marginalize over the calibration and beam uncer- 
tainty in order to account for remaining systematic uncertainties 
in our results. To provide an extra large-angular-scale constraint 
we combine our data with the results from the COBE-DMR 
satellite as provided by Gorski et al. (1996). 

Within the chosen family of inflationary models we put strin- 
gent 95% confidence level constraints on the total density, £1 = 
0.9^'jg, baryon density, VL\,h 2 = 0.033 ±0.013 and power spec- 
trum normalization Cio = 690^j25^ 2 . Our 95% confidence 
limit on the cold dark matter density is Q c d m h 2 = 0.17^'^. 
However, this result is mostly determined by the priors defin- 
ing the database parameter range, as discussed in Jaffe et al. 
(2001). We find a strong degeneracy between the optical depth 
to the last scattering surface, r c , and the primordial power spec- 
trum index, n s . In this case, the degeneracy restricts the param- 
eters to a subspace allowing us to derive a combined constraint: 
n s ~ (0.99±0.14) + 0.46t c (95% c.l.), for r c < 0.5. We note 
that we recover the 95% upper limit ( 0.4) on r c , derived by 
Griffiths, Barbosa, & Liddle (1999), if we constrain the spectral 
index to be < 1.2 as was assumed by those authors. Indepen- 
dent of the value of the optical depth, we can put a very firm 
^ 99% lower limit on the spectral index, n s > 0.8. Assuming 
no reionization (t c = 0), we can get both lower and upper limits 
on the spectral index reading « s = 0.99 ±0.14. Alternatively, 
fixing n s at unity gives us a 95% c.l. upper limit r c < 0.26. 

The x 2 °f the best fit model is 30 for all 41 points 
used in the fittings, and 4 for the 13 points of MAXIMA-I 
only. The best fit model parameters (ilb,^cdm)^A,T c ,M s ,/z) = 
(0.07, 0.68, 0. 1 , 0.0, 1 .025 , 0.63) are characterized by high mat- 
ter and low vacuum energy content (see also Balbi et al. 2000). 
However, due to strong degeneracy between il m and fl\ (e.g., 
Zaldarriaga, Spergel & Seljak 1998) we can easily find models 
comfortably fulfilling both MAXIMA-I+DMR and supernovae 
constraints (see Figure 3). These data sets constraint jointly 
the matter and vacuum energy densities to be: (Q m ,flA) = 
(0.32! 1 1 4,0.65! ;^)(95%c.l.). 

The constraint on il b h 2 mentioned above is compatible 
with the best determination to date of the baryon density 
based on measurements of primordial deuterium and calcu- 
lations of standard BBN (Buries, Nollett, & Turner, 2001, 
Tytler et al. 2000), n b h 2 = 0.020 ± 0.002. The consis- 
tency of the data with BBN becomes even more apparent 
by constraining our parameter estimation on the BBN value 
of ttbh 2 , which we approximate by fixing Sl\,h 2 at the grid 

ses of freedom; the quotoed numbers are very close to the maxima of those 
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value nearest to the BBN prediction (i.e., Q\,h 2 = 0.02). The 
best fit model then has parameters {il^,n c ^ m ,fl\,T c ,n s ,h) = 
(0.07, 0.78, 0.0, 0.0, 1 .0, 0.53) with X 2 - 1 for the 13 MAXIMA-I 
points only. This very good fit emphasizes the compatibility of 
our data with other cosmological measurements. The spectrum 
of this model is indicated in Figure 1 with a dotted line, and it 
shows that even in that case, the best fit model has a higher third 
than second peak, if £(£ + 1)Q is plotted versus I. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Fig. 3. — Constraints in the f2 ra -f2A plane from the combined 
MAXIMA-I and COBE-DMR datasets. The shown contours corre- 
spond to 68, 95, & 99% likelihood-ratio confidence levels. The 
bounds obtained from high redshift supernovae data (Perlmutter et al. 
1999; Riess et al. 1998) are also overlaid as well as the confidence 
levels of the joint likelihood. 

We also compute constraints on the ratio of bin powers, TZo\, 
as imposed by the MAXIMA- 1 data within the discussed family 
of models. We find that the most likely value of that ratio is 
0.68 j'jZ at the 95% confidence level, consistent with (though 
on average higher than) our result in Sect. 2 (cf., Figure 2). 

4. IMPLICATIONS FOR COSMOLOGICAL MODELS 

Inflation-based models provide us with an abundance of ex- 
cellent fits to the MAXIMA-I extended power spectrum. These 
include the best fit models found by the analysis of the first data 
sets of BOOMERANG-LDB and MAXIMA-I (Jaffe et al. 2001, 



Balbi et al. 2000). Constraints on cosmological parameters us- 
ing DMR and MAXIMA-I confirm the near flatness of the Uni- 
verse and, when combined with recent supernova data (Riess et 
al. 1998, Perlmutter et al. 1998), support the need for the non- 
zero vacuum energy density. The most likely baryon density 
indicated by our data is found to be somewhat higher than the 
preferred BBN value, yet the latter is within the 95% confidence 
range of our determination. Moreover, there are excellent fits 
to the data with the baryon fraction at the BBN value. 

Our results constrain the power at £ ~ 860 to be ^ 1 700/iA' 2 
on the 95% confidence level. They also indicate an increase in 
the power spectrum between £ ~ 410-785 and £ ~ 786-925, 
yielding the ratio of the corresponding bin powers equal to 
75-01 = 0-49_ 2 1 in a general case, or 1Zq\ = 0.68^'^ within 
the considered family of inflationary models. We have shown 
that such constraints can be easily fulfilled by these inflationary 
models. However, they impose strong requirements on some 
non-standard models which were found to better accommodate 
the first results of the BOOMERANG-LDB and MAXIMA-I ex- 
periments. Those models, by design, have an amplitude of 
power in the intermediate range of £ below the typical predic- 
tions of the standard inflationary models. However, they also 
tend to have lower power at the high £ end and therefore, to 
match our new constraint at £ ~ 860, they would require even 
higher baryon abundance than the family of the models con- 
sidered in Sect. 3. These disfavored models include: the mix- 
ture of inflation and topological defect models (Bouchet et al. 
2000, Contaldi 2000), hybrid adiabatic plus isocurvature mod- 
els (Stompor, Banday, & Gorski, 1996, Enqvist, et al. 2000), 
broken power-law primordial spectra models (Griffiths, Silk, & 
Zaroubi, 2001, Barriga et al. 2001), or models with a delayed 
recombination (Peebles et al. 2001). Further investigation is 
required to fully elucidate status of these models. 
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